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ABSTRACT 

One of the first stages of planet formation is the growth of small planetesimals and their accu- 
mulation into large planetesimals and planetary embryos. This early stage occurs much before the 
dispersal of most of the gas from the protoplanetary disk. Due to their different aerodynamic prop- 
erties, planetesimals of different sizes and shapes experience different drag forces from the gas during 
this time. Such differential forces produce a wind-shearing (WISH) effect between close by, different 
size planetesimals. For any two planetesimals, a WISH radius can be considered, at which the differ- 
ential acceleration due to the wind becomes greater than the mutual gravitational pull between the 
planetesimals. We find that the WISH radius could be much smaller than the gravitational shearing 
radius by the star (the Hill radius). In other words, during the gas-phase of the disk, WISH could 
play a more important role than tidal perturbations by the star. Here we study the WISH radii for 
planetesimal pairs of different sizes and compare the effects of wind and gravitational shearing (drag 
force vs. gravitational tidal force). We then discuss the role of WISH for the stability and survival of 
binary planetesimals. Binaries are sheared apart by the wind if they are wider than their WISH radius. 
WISH-stable binaries can also inspiral, and possibly coalesce, due to gas drag. Here, we calculate the 
WISH radius and the gas drag-induced merger timescale, providing stability and survival criteria for 
gas-embedded binary planetesimals. Our results suggest that even WISH-stable binaries may merge 
in times shorter than the lifetime of the gaseous disk. This may constrain currently observed binary 
planetesimals to have formed far from the star or at a late stage after the dispersal of most of the disk 
gas. We note that the WISH radius may also be important for other processes such as planetesimal 
erosion and planetesimal encounters and collisions in a gaseous environment. 
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1. INTRODUCTION 

Gravitational encounters between planetesimals play 
an important role in the evolution of proto planetary disks 
and planet formatio n (e.g., reviews by iLissauei) [19931: 
iGoldreich et al.ir2004D . Planetesimal growth likely occurs 
while the planetesimals are still embedded in a gaseous 
disk. Studies of gas-planetesimal interactions have shown 
that gas can affect the velocity dispersion of planetesi- 
mals (e.g. .Nelson fc Gressel,2010 ). may help the forma- 
tion of large planetesimals thro ugh clumping of planetes- 
imals ([Chiang fc Youdiiil 120101 and references therein), 
and can lead to fa st inspiral of planetesimals into the star 
throu gh gas-drag ([Nakagawa et al.iri986l : IWeidenschillind 
119771 and references therein). Here we focus on a differ- 
ent aspect of planetesimals embedded in a gaseous disk, 
namely the close interaction between pairs of single plan- 
etesimals in a gaseous environment. 

Planetesimals likely vary in size and shape, and there- 
fore have a wide range of aerodynamical properties, 
which affect their interaction with surrounding gas. In 
particular, planetesimals of different sizes and/or shapes 
experience different drag forces from the head wind 
they encounter in the gaseous disk. The difference be- 
tween thB_Jorces_a£tiiig_on two different-size planetesi- 
mals (| Weidenschillinei[l 9 771) can change their relative tra- 
jectories with respect to the ir unperturbed motion in the 
absence of gas (for example lOrmel fc Klahiir2010l consid- 
ered planetesimal interactions in gas rich environment; 
their study focused on planar encounters and drag law 
regimes which are linearly dependent on velocity). 

During an encounter between two different-size plan- 



etesimals, the different forces experienced by the two 
components as a result of gas drag generate a wind- 
shearing (WISH) effect, which could be stronger than 
their gravitational interaction. For any two planetesi- 
mals, we consider the radius, which we term the WISH 
radius, at which the differential acceleration due to aero- 
dynamical wind-shearing becomes greater than the mu- 
tual gravitational pull between them. In the following we 
explore this new distance scale and discuss its implica- 
tions, including the stability of binary planetesimals. In 
addition, we study the evolution of WISH-stable binary 
planetesimals in gas. Such binaries dissipate their orbital 
energy through gas drag and may inspiral to form closer 
binaries or even coalesce during the typical lifetime of a 
protoplanetary disk. 

We begin by deriving the WISH radius and discussing 
the effects of gas drag on particles of different sizes. We 
then calculate the WISH radius for two planetesimals of 
arbitrary effective sizes (Section[2|). In Section[31 we con- 
sider the evolution of binary planetesimals embedded in a 
gas disk, including WISH stability (Section [3Tt and gas 
drag-induced inspiral (Section 13. 2p . Finally, we discuss 
various other possible implications of our results (Section 
[4]) and summarize (Section [5]). 
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2. GAS DRAG AND THE WIND-SHEARING RADIUS IN 
PROTOPLANETARY DISKS 

Planetesimals of different sizes embedded in the same 
gaseous environment experience different drag forces and 
hence different accelerations. The differential accelera- 
tion between two planetesimals of mass rrib and TOs due 
to the wind-shearing effect is given by 



rrib 



3pp 
47r 



Foirb) Fairs) 



(1) 

where Fo is the force exerted on a particle due to gas 
drag. Throughout this paper, we perform our calcula- 
tions for spherical particles of constant density, pp, so 
that the mass of a planetesimal with radius r is m = 
(4/3)7rppr^. In Equation ([T]), planetesimal masses nib 
and nis correspond to radii rb and r^, respectively. Real 
planetesimals could have different aerodynamical prop- 
erties (e.g. they may not be spherical and/or they could 
be porous); however, calculations analogous to those pre- 
sented here may be performed for any form of Fjj (m) . 

For small separations over which the environmental 
conditions (gas density and temperature) are approxi- 
mately the same, the differential WISH acceleration be- 
tween any two planetesimals is independent of the dis- 
tance between them. In this case, Aaws can be used 
to define an important distance scale, which we term 
the WISH radius. To define this scale, we adopt a simi- 
lar approach to that used to define the the gravitational 
tidal-shearing radius, i.e., the Hill radius. 

The Hill radius (sphere) is the distance at which the 
gravitational influence of a planetesimal or a planet with 
mass m and radius r, orbiting a star with mass AI^, at 
radial distance a, becomes comparable to the tidal per- 
turbation by the star. It is given by 
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where the second expression is derived for a spherical 
planetesimal with m ^ M^. A test particle located close 
to the Hill radius, or beyond it, is strongly affected by the 
gravitational pull of the star. If it begins in orbit around 
the planetesimal, its orbit will be perturbed and is likely 
to become unstable. The exact distance up to which a 
binary orbit can remain stable also depends on its orbit 
direction, e.g. prograde or retrograde with respect to the 
orbit of the planet around th e star plamilton fc BurnsI 
IT99lt IShen fc TremamilMl IPerets fc Naozll2009D . or 
more generally, the relative inclination of the particle's 
orbit. In the following we adopt the simple deflnition 
given by Eq. ([2]) . 

Following the definition of the Hill radius we can now 
define the WISH radius. This radius is defined as the 
distance between two planetesimals for which the differ- 
ential WISH acceleration between them equals their mu- 
tual gravitational pull. Equating /S.aws with the gravi- 
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tational acceleration a 

separation dbin between the two planetesimals equal to 
the WISH radius, which we define as 
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i Gjmb + nis) 
Aaws 
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Beyond this limiting radius even two planetesimals which 
are formally gravitationally bound (in the absence of 
WISH) would be sheared apart by the wind. 

In order to calculate the specific value of the WISH ra- 
dius for any given pair of planetesimals, we first need to 
understand the gas-drag force applied on planetesimals 
which face a head wind. This depends on the specific 
regime of the gas-planetesimal interaction, since differ- 
ent gas-drag laws apply under different conditions. We 
review gas drag laws in Section [2. 1) calculate Aaws ex- 
plicitly for several regimes in Section 12. 2| and combine 
these to provide self-consistent calculations of Rws as a 
function of planetesimal size in a fiducial disk (Section 
[231) . 

2.1. Drag laws 

The appropriate gas-drag force on a planetesimal of 
radius r moving through gas at relative velocity Vrei de- 
pends on r/A, where A = fj,/{pga) is the mean free path 
of the gas, a is the cross-section for gas-gas coUisions, 
and fj, is the mean molecular weight. For planetesimals 
with diameters larger than the mean free path of the 
gas, the drag force also depends on the fluid Reynolds 
number Re — 2rvrei/(0.bvth^)- Here, 2r is the diame- 
ter of the planetesimal. The gas has kinematic viscosity 
{l/2)vth^i temperature T, and mean thermal velocity 
(for a Maxwelhan distribution) vth — (8/7r)^/^Cs, where 
Cs = (kT/ ^^/"^ is the sound speed and k is Boltzmann's 
constant. The various gas drag r egimes can be summa- 
rized as follows (where we follow iWeidenschillingl II977L 
who in turn follows I Whi^HIi ll 9731 ) Fl 

When r < A, drag may be modeled by considering 
individual and independent particle collisions, and the 
Epstein regime applies (for subsonic Vrei, which is ap- 
propriate for our problem): 



Fd = -^T^PgVthVreir 



(4) 



For r > A, the gas must be modeled as a fluid. We take 
r = (9/4)A as the boundary between these regimes. At 
low Reynolds number, the gas/particle boundary layer 
dominates (Stokes drag), while at high Re, the gas exerts 
a Ram pressure force on the particle, so that the drag law 
is 



Fo=3TrpgVthVrei>^r 



for i?e < 1 
for Re > 800 



Stokes 
Ram. (5) 



An intermediate regime exists for 1 < Re < 800. A/lore 
generally, the full range of Reynolds numbers can be fit- 
ted with a drag law of 
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where CoiRe) can be fitted with an empirical formula 
based on r ecent experimental data in the regim e 10^^ < 
Re < 10^ (H rown fc Lawleii[200l lChend[2009t compati- 

^ Planetesimals move through the protoplanetary disk at sub- 
sonic velocities. For Vj.^i > Cg, ram pressure drag applies (c.f. 
Equation [5j . 
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ble with older data used bv lWhipplill973[ ). yielding 
24 
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(l + 0.27i?e) 



0.43 



0.47[1 - ea;p(-0.04i?e" -^**)] 



(7) 

We use Equations ^ and ([7|) for our drag law in the 
calculations that follow, except when a single, specific 
drag law is specified, in which case we use Equations (|4]) 
and dni). 

2.2. The wind- shearing differential acceleration 

When two particles with radii rf, and r^ experience gas 
drag in the same drag regime, we may obtain a simple 
expression for the wind-shearing differential acceleration, 
/S.aws- For example, in the Epstein regime {rs,rb < A), 



rlvrei{rb) r1vrei{rs) 



rub 



Epstein. 
(8) 

In general, the relative velocity between each of the plan- 
etesimals and the gas could differ, in which case even 
planetesimals of the same size can experience a differen- 
tial WISH acceleration. For bound binary planetesimals 
(which we discuss in Section 13.11) the velocity relative to 
the gas of the two components should be, on average, ap- 
proximately the same, so that Vrel = Vrelifb) = Vrelifs)- 

In this case, Eq. (jH]) simplifies further: 
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where the last expression is for rt ^ r^. Similarly, under 
the same assumptions, the Stokes regime produces 



AawS - iT^PgVthVrel>^ 



P VthVrel 



4 ppcr rj 
and in the Ram pressure regime. 
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We note that in the Stokes regime, Aaws does not de- 
pend on the density of the gas. 

More generally, the differential acceleration can be ob- 
tained accurately for any combination of two planetesi- 
mals in different drag law regimes and moving through 
the gas at different velocities. This can be done by using 
Eqs. (gl) and ^ or Eqs. ^ and (O to calculate the 
appropriate drag force on each planetesimal. 

2.3. The wind-shearing radius 

Given the expressions in Equations ([9))- (ITl1) . we can 
obtain relatively simple formulas for the wind-shearing 



radius (Equation [3]) for two planetesimals in the same 
drag law regime with Vrei{i'b) — Vreiir^) and rt, 3> 
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Epstein, 
Stokes, 
Ram, 



in the Epstein, Stokes, and Ram pressure regimes, re- 
spectively. In fact, these expressions apply as long as the 
gas accelerates the smaller body more effectively than the 
larger body, even if their drag regimes are different. In 
this case the drag regime in Eq. would correspond 
to that of the smaller body. 

More generally, we can calculate Rws for two plan- 
etesimals of arbitrary size (in the same or in different gas- 
drag regimes) as a function of the properties of the gas 
in which they are embedded. When < Rws < Rh, 
the WISH radius represents the limiting separation of a 
binary planetesimal (or a pair of small satellites) in a 
gaseous environment. 

In Figure [U we show the calculated Rws and the re- 
sulting binary stability radius for a planetesimal with ra- 
dius rb = 10 km, orbited by a smaller body with a range 
of sizes. This calculation is performed at 1 AU from a 
solar-mass star in a disk having the following parameters. 
We choose a disk temper ature of T = ro(a/AU)~"^/^ with 
To = 120 K, follo wing iChiang fc Youdin (2010), who 
adapt the results of lChiang Sz Goldreich (1997.) for a disk 
around the young Sun. Varying the value of Tq within a 
reasonable range for Sun-like stars does not qualitatively 
change our results. We take the surface density of the 
disk to be Eg = So(a/AU)"\ with Eq = 2 x 10^ g/cm^. 
This choice is roughly consistent with the minimum-mass 
solar nebula at 1 AU as well as with observed dust sur- 
face density profiles at distances larger than ~ 20 AU 
in extrasolar disks, taking a dust to gas mass ratio of 
1 : 100 (typically assumed in the modeling of protoplan- 
etary disks, e.g. Andrews et al. 2010). Protoplanetary 
disks likely exhibit a range of surface density profiles 
across different systems and at different times within the 
same system. We discuss the impact of varying the gas 
surface density in the Appendix. 

Given the above choices, the disk scale height H is 
given by 



- - ^ - 0.022 (^) 
a fla \AUJ 



2/7 



(13) 



where fl — {GMi,/ a'^Y^'^ is the Keplerian orbital fre- 
quency. The gas volume density profile is then 



-16/7 



l/cm. (14) 
Using the neutral collision cross-section a ~ (SA)^ ~ 
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Figure 1. The binary stability radius for a r;, = 10 km planetes- 
imal with a companion of radius rs < r^, at 1 AU from the star 
in our fiducial disk. At radii between the physical size and the 
Hill radius of the large planetesimal (lower and upper dashed lines, 
respectively), the stability radius equals the WISH radius, Rws 
(see text). The planetesimals are moving at a relative velocity of 
w^j,; ~ 0.5c^/i)jf with respect to the gas. The relevant drag regime 
for the small body changes with Ts , and the WISH radius is well 
approximated by Equation I I12II . A pair of bound planetesimals (a 
binary) can only exist in the shaded region; it can not reside lower 
than the physical size (collision), above the Hill radius (gravita- 
tionally unbound by tidal shearing from the star) or to the left of 
the WISH radius (sheared apart by the wind). 



10 cm^, the mean free path of the gas is 

A cm 

Uga \AUJ 



(15) 



where the gas number density Ug — Pg/ii and we have 
used fj, — 2.3mH with mn equal to the mass of a hydro- 
gen atom. 

For these calculations, we assume that the relative ve- 
locity between the binary and the gas is equal to the ve- 
locity of a single planetesimal with radius equal to that 
of the larger component of the binary, r^, as it moves 
through a smooth disk under the influence of gas drag. 
This approximation is valid for 3> rs- In reality, a 
bound binary will move through the gas at a velocity 
that r eflects th e drag on bot h binary components. Fol- 
lowing lYoudinl (2010; see also iNakagawa et all ll 9861 ). we 
set the relative velocity between a planetesimal and the 



gas to be v. 



rel 



(v'^ , , ^)^/^ with 
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(16) 
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with T] = [vk — Vg,4>)/vK, so that rjVK equals the differ- 
ence between the azimuthal gas velocity, Wg.^, and the 
Keplerian velocity, vk = ^ol- We use the approximate 
value rj ~ 0.5(Cs/f|-). We calculate the stopping time 
tg = mVrei/Fu and relative velocity v^ei of a planetes- 
imal iteratively, using the drag law in Equation ([71), in 
order to achieve self-consistent values for these and hence 
for Fr) in all drag regimes. Note, however, that our choice 
of Vrei represents the velocity of a planetesimal moving 
through a uniform disk and does not take into account 
turbulence. In addition, even in a smooth disk planetes- 



imal growth likely occurs in regions of enhanced solids, 
which may accelerate the disk gas to more nearly Kep- 
lerian speeds, reducing this relative velocity. We discuss 
how our results vary as a function of relative velocity in 
the Appendix. 

Under our assumed conditions, a 10 km planetesimal 
at 1 AU orbits at approximately the Keplerian velocity, 
so that in Figure [TJ Vrei ~ 0.5cI/vk- At this relative 
velocity, the Reynolds number for a planetesimal with 
radius r, is 



Re=- 



2r,ti, 



2 



0.5ut/iA 

0.02 f^) f-^ 

VI cm/ \AU 



(18) 



The transition from the Epstein to the Stokes regime for 
the small planetesimals may be clearly seen in Figure [1] 
as a change in the slope of the WISH radius from 1/2 to 
unity. This behavior is matched by Eq. (fT2|) . Though 
the large body in this plot always remains in the Ram 
pressure drag regime, it is not accelerated much by the 
gas, and Aaws is dominated by the acceleration of the 
small companion. The agreement between these results 
and Eq. (fT^ reflects the fact that in this regime, the 
small companions are accelerated more effectively by the 
gas than the 10 km large body. 

Figure[2]displays the calculated WISH radius as a func- 
tion of the small planetesimal size for various sizes of the 
large planetesimal and at different distances from the 
star. Also shown for comparison are the physical size of 
the big planetesimal and its Hill radius. The WISH ra- 
dius diverges for equal size planetesimals (with the same 
velocities relative to the gas), since they experience the 
same gas drag, and their differential WISH acceleration 
approaches zero. 

3. BINARY PLANETESIMALS IN A GASEOUS 
ENVIRONMENT 

A non-negligible fraction of currently observed plan- 
etesimals in the Solar system (including asteroids and 
Trans-Nept unian objects) are found to be members of b i- 
naries (e.g.. lRichardson fc Walshll200fit l?foirern] [2008^ 
Binary planetesimals can teach u s about the dynami 
cal evolution of the Solar system (iPerets fc N aoz 200_^ 
Murrav-Clav fc Schhchtingj 120111: iParker fc'Kavelaar:,. 

and can play a role in planet formation an d 
planetesimal growth (jNesvornv et al.ir2010t lPeretsl[2010( ) . 
Study of the interactions of binary planetesimals with gas 
is therefore important for understanding the formation, 
stability and evolution of these binaries and their impli- 
cations. In the following we discuss the effect of WISH 
and gas drag inspiral and coalescence of binary planetes- 
imals in gas. 

3.1. Wind-shearing disruption of binary planetesimals 

We have already alluded to an immediate consequence 
of the wind-shearing radius for binary planetesimals, 
namely that it provides a new stability criterion for their 
survival (see Figure [1]). In a gas free environment, bi- 
nary planetesimals are stable as long as their separation 
is smaller than the Hill radius, whereas wider binaries 
are destabilized and disrupted by the tidal gravitational 
shearing from the star. However, in the presence of gas, 
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Figure 2. The WISH radius, Rws, (see text and Figure [T]l for 
two planetesimals of sizes rs and ri, as a function of , at dis- 
tances of 1, 5 and 40 AU from the star (top, middle, and lower 
panels, respectively). The plotted lines show Rws foi' planetesi- 
mals with radii ri, = 10~^ — 10* cm, in logarithmic jumps. The 
shaded regimes show the regions where the WISH radius is smaller 
than the physical radius or larger than the Hill radius of the big 
planetesimal. We do not show lines corresponding to big planetesi- 
mals for which the WISH radius is smaller than the physical radius 
unless the smaller planetesimal nearly equals the larger planetesi- 
mal in size (e.g. r;, < 1 m at 1 and 5 AU). 

the Hill radius stability limit should be replaced by the 
WISH radius when Rws < Rh (binaries wider than the 
Hill radius are always unstable). We find that the sta- 
bility criterion for binaries embedded in gas is 



dbin < min{RH, Rws)- 



(19) 



Because collisions prevent binaries from forming with 
dbin < "Tb, no stable binaries are possible when Rws < i"b- 
Given these considerations, the limiting separations of 
binary planetesimals in our fiducial disk as a function of 
size and distance from the star may be read from Fig- 



ure [5] (see also Figure [T]). For planetesimal sizes spanning 
a wide range, < Rws < Rh in our fiducial disk and 
this limiting separation is equal to the WISH radius. Bi- 
nary planetesimals can therefore be strongly affected by 
WISH, most notably for smaller planetesimals closest to 
the star. Generally we find that for a wide range of bi- 
nary and disk properties the WISH radius determines 
the stability rather than the Hill radius. A gaseous envi- 
ronment qualitatively changes the spatial dependence of 
binary stability in a protoplanetary disk, as the spatial 
dependence (distance from the star) of the Hill radius 
and that of the WISH radius differ. This can be seen 
in Figure [31 which compares the WISH radius with the 
Hill radius for a given binary pair as a function of the 
distance from the star (also compare the panels in Fig- 
ure [2]). We note that binaries with radii of ~ 100 km 
and components of roughly equal mass, comparable to 
many observed Trans-Neptunian and asteroid binaries, 
have their stability determined by the Hill radius at all 
distances from the star in our fiducial disk. 

More details regarding the dependence of our results on 
the gas density and the relative velocity of planetesimals 
with respect to the gas can be found in the Appendix. 

3.2. Gas drag-induced inspiral of binary planetesimals 

Many studies have demonstrated that single planetesi- 
mals inspiral toward their host star du e to gas drag (e.g., 
Weidcnschilling 1971; for a review see lChiang &: YoudinI 
20101 and references therein) . A similar process can cause 
the inspiral of a binary planetesimal into a closer mu- 
tual orbit, possibly ultimately leading to coalescence. In 
the following we explore the evolution of the mutual or- 
bits of binary planetesimals in gas, and we provide the 
timescales for their coalescence. 

For simplicity we restrict our discussion to binary plan- 
etesimals with rris rrib, but our results can be simply 
generalized to arbitrary mass ratios. We study the evolu- 
tion of a binary orbiting the Sun in the plane of the pro- 
toplanetary disk. For simplicity, we assume throughout 
that the mutual binary orbit begins circular and dissi- 
pates orbital energy on a timescale much longer than the 
mutual orbital period. The two components of the binary 
therefore orbit one another on roughly circular trajecto- 
ries at a (shrinking) binary separation of dbin- In some 
circumstances binaries can evolve on faster timescales; 
such short term evolution relates to binary planetesimal 
formation through gas dissipation and may result in fast 
coagulation of planetesimals. These latter processes will 
be discussed elsewhere (see Murray-Clay & Perets, in 
preparation). 

In this limit, the binary loses angular momentum L 
on a timescale of \L/L\ — msVbin/ {Fd), where Vbin 
is the binary orbital velocity of the small body and 
(Fd) is the gas drag force on the small body, aver- 
aged over a binary orbital period. Equivalently, the bi- 
nary loses orbital energy E on a timescale of \E/E\ — 
O.bmsvl^^/ {{F]:)) Vbin)- Inspiral therefore proceeds on a 
timescale Tniorgc = dbin/dbin = (1/2)-^/-^ = E/E, so that 



dbin 
dbin 



1 msVbi, 

2 (Fd) 



(20) 



In calculating the gas drag force, we must average over 
an orbital period because the center of mass of the binary 
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Figure 3. Comparison of the WISH radius and the Hill radius 
as a function of planetesimal size and distance from the star. The 
WISH radius restricts binary stability most at larger separations, 
where the tidal force from the star becomes weak. Top panel: 
The planetesimal size for which Rws('''b) = Rnifb) (dotted line), 
Rws = 0.1-Rff (dashed line) and -Rvi's(^b) = (solid line), as 
a function of the distance from the star, for a binary planetesimal 
with size ratio rg/r^, = 0.1. Middle panel: Same a top panel, but for 
^'s/^'h = 0.5. Larger size ratios reduce the impact of WISH, making 
binaries more stable. Binaries in the bottom shaded region can not 
survive WISH and will be disrupted by the wind. The stability of 
binaries in the top region (above the Rw = Rh line) is set by the 
Hill radius rather than by the WISH radius, i.e. in that region 
gravitational shearing is stronger than WISH. Bottom panel: The 
binary WISH stability radius as a function of the distance from the 
star for three binaries with a size ratio Ts/r^, = 0.5, and primaries 
of = 1 km (top red solid line), r^, = 0.1 km (middle black 
solid line) and rj, = 0.01 km (bottom blue solid line). Also shown 
(dashed line) is the spatial linear dependence of the Hill radius, 
corresponding to the stability radius for a non-gaseous planetesimal 
disk. Note that Rul^h is independent of r^. 



is m.oving with respect to the background gas with rela- 
tive velocity v^isk as it orbits the Sun. Given rtis ^ m^, 
the relative velocity between the binary and the disk gas 
is approximately the velocity at which the large body 
would move through the gas on its own, given by Equa- 
tions (Uni) and (fT7|) . 

We now provide analytic expressions for the infall time 
in two different gas-drag regimes, the regime which is lin- 
ear in velocity (corresponding to the Stokes and Epstein 
regimes; Section [3. 2 .ip and the quadratic (ram pressure) 
regime (Section l3.2.2p . In the quadratic regime, the type 
of evolution depends on the ratio Vbin/vdisk- In practice 
more complicated regimes exist (see Section l2.ip . which 
we integrate numerically for parameters relevant to plan- 
etesimals in a protoplanetary disk in Section 13.2.31 



3.2.1. Linear 



regime 



In the following treatment, we assume that v^in re- 
mains constant over a single binary orbital period Pbin, 
which is good for Vbin/vbin ^ ^bin- Note that this as- 
sumption requires not only that Tmerge S> fbin/2 but also 
that msVbin/FD,disk » -Pbin, where FD,disk is the drag 
force experienced by the small body moving at relative 
velocity Vdisk with respect to the gas. We address the 
complication of non-circular orbits in future work. 

In the linear regime, Fd cx VreU with v^ei equal to the 
relative velocity of the small body with respect to the 
gas, containing components from the binary orbit and 
from the overall motion of the binary through the gas 
disk. Therefore Fo,i = Fo/vrei is constant over the 
binary orbit. The linear regime is valid for the Epstein 
and Stokes drag regimes, but the value of in the two 
regimes differs (see Sect ion [^?T|) . We may now express the 
orbit-averaged drag force as 



{Fl 



1 

2^ JO 
Fds 
27r 



Fode 

2tt 



(21) 



{vbin sin^ + Vdisk)d9 = Fo^iVbi, 







where 9 is the angle of the binary in its orbit. The 
term Vbin sin 9 is the bulk velocity component of the small 
planetesimal parallel to the the direction of motion in the 
binary frame of reference, so that v^ei = Vbm sin 9 + Vdisk ■ 
Over a full orbit the contribution from Vdisk averages out 
and 



t 



stop 



(22) 



with tstop equal to the stopping time of a single small 
planetesimal in the gaseous protoplanetary disk: 



istop — 



Pp 



PgJ Vth 



Fd; 



4 /pp 

9 \P9 



Epstein 



Stokes. 



'th 



Recall that in the linear regime, the stopping time is inde- 
pendent of the relative velocity between the planetesimal 
and the gas. Note that single planetesimals with stop- 
ping times longer than an orbital time inspiral into the 
star on a timescale of '^tstop/v- The same processes are 
at work in both cases — infall into the star is slower than 
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binary coalescence because the gas and planetesimals or- 
bit the star together, reducing their relative velocities. 

The timescale for coalescence is independent of dbim 
and the total merger time for a binary is 



J. rr 



1 I'^o 
-In I — 



(23) 



where dbm = do initially, and rf, is the final binary sepa- 
ration before coalescence. 

3.2.2. Quadratic (ram pressure) regime 

We now consider the quadratic regime, for which cx 
"f^re/' appropriate for ram pressure drag. Following the 
same procedure as above, but using Fd_2 = Fd/vI^j with 
Fd,2 a constant, we get 



D 



Fd^2 

27r 



2tt 



iVbin sm ( 



vdiskTde 



-FD.2Vbin 



1 

1+2 



Vdisk 



(24) 



In other words, the ram pressure drag force re- 
quires an effective relative velocity correction of [1 -I- 
0-5{vdisk/vbin)'^] — in this case the contribution from the 
bulk velocity drag did not average out. 
Now, 



1 + 0.5{vdtBk/vbiny 



(25) 



where tstop(vbin) is the stopping time for Vrei = Vbin- In 
the quadratic regime, tstop is not independent of Vrei, so 
to make dependences clearer, we rewrite this expression 
as 



D.2) 



Vbin[i + O-H^dtsk/Vbin)^] 

rUs 

TTT^ , Vbin ^ Vdisk 

^-rD,2Vbin 

PoWdtsk 



(26) 



, Vbin ^ Vdisk 



Plugging in Fd^2 for ram pressure drag and Vbv. 
{Guib/dbinY^'^, this corresponds to 



Pp_ 

0.66 V Pg 



'^ulNGrtib 



, Vbin ^ Vdisk 



. '^y/Gmb/{dl'^V%^^) , Vb;n < Vd^sk 

Integrating, we find a total merger time of 



(27) 



Pp 



merge ~ r> oo \ 

0.33 \Pg, 



X < 



y/GlTLb 



, Vbin ^ Vdisk 



2VQrrib 



(28) 



disk 



,1/2 



, Vbin ^ Vdisk 



The merger time in the ram pressure regime depends 
strongly on the ratio between the binary mutual orbital 
velocity and its bulk velocity around the star. When 
Vbin ^ Vdisk, the merger proceeds most slowly when dbin 
is largest, while for Vbin ^ Vdisk, the final coalescence at 
dbin Tb takes the longest time. 

3.2.3. Implications 

As can be seen in our analytic derivation, the timescale 
for the inspiral of a binary planetesimal is dependent on 
its environment and on the binary properties. Figure |4] 
shows the binary merger timescale as a function of small 
planetesimal size for a range of big planetesimal sizes and 
separations from the star. To make this figure, we cal- 
culate the integral {Fd) — J^^ FD{vrei)d9 numerically 
with Vrei = Vbin sin + Vdisk , using the full expression 
for Fd embodied in Equation ([7]). This, for example, 
allows the relevant drag law to vary as a function of 9 
if appropriate. We choose either the merger timescale 
evaluated at dbin = min{RH,Rws) or at dbin = fb, 
whichever is larger. We maintain our assumption of cir- 
cular orbits. This assumption is only valid for merger 
timescales longer than of order the orbital period of 
the binary around the Sun, Porb- The assumption that 
^^sVbin /FD,disk ^ fbin breaks down for binaries with 
large bodies smaller than ~100m-lkm in size, making 
our (already short) merger timescales upper limits in 
these cases. 

We find that binary planetesimals over a wide range of 
masses inspiral and likely merge in times much shorter 
than the typical lifetimes of gaseous protoplanetary 
disks. Binary asteroids with components having radii less 
than a fe w tens of km have been ob served in the main 
belt fe.g.. iRichardson fc Walshll2006l ). Such pairs were 
not likely to survive for long in a gaseous disk (see top 
and middle panels of Figure HJ. Observed binary TNOs, 
however, have radii larger than a few tens of km (e.g., 
Noll et al. 2008), and could have survived for more than 
a Myr (bottom panel of Figure [4}. If binary minor plan- 
ets formedjii_the_prm gas disk (e.g. as suggested 
bv lNesvornv et ani2010l ). those with small components 
are less likely to have survived to this day. Given this 
formation scenario, the orbital characteristics of even the 
largest binary asteroids likely changed due to their evolu- 
tion in gas. The currently observe d orbital properties of 
binary asteroids (|Naoz et al.ll2010t) are therefore unlikely 
to reflect only their properties at birth; e.g. binaries 
born with wide separations were likely to inspiral into 
more compact configurations. While the orbits of cur- 
rently observed binary TNOs were likely unaffected, as 
binary TNOs with smaller components are found, this 
effect will need to be considered. 

4. OTHER ASPECTS OF GAS-PLANETESIMAL 
INTERACTIONS 

As we suggested above, the WISH radius can have im- 
portant implications for gas-planetcsimal interactions. In 
this study wc mainly focused on the implications of the 
WISH radius for binary planetesimals. However, this 
scale could be important for other processes, similar to 
the role played by the Hill radius in gas-free environ- 
ments. Here we only briefly mention these issues, which 
will be discussed in detail (and more quantitatively) else- 
where. 
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Figure 4. Merger timescales, Tmerge for binaries located 1 (top), 
5 (middle), and 40 (bottom) AU from their host stars, as a function 
of the radius of the smaller binary component, r^. Binary plan- 
etesimals with a wide range of masses inspiral and likely merge 
on timescales shorter than the typical lifetimes of protoplanetary 
disks (gray region). Colors correspond to different large planetes- 
imal sizes, matching the colors in Figure [5] For large planetesi- 
mals with radii < lO'^km, the merger time is independent of 
ri,. The merger time is reduced for rj, = lO^km (orange, lower 
curve) as the ram pressure drag regime becomes important. For 
ri, < 100m (top) or 1km (middle and bottom), our assumption 
that rUsViyin/ F]j iHgk S> fbin breaks down, and our curves are up- 
per limits. Timescales shorter than the orbital period around the 
star, Porbi (hashed region) indicate that inspiral will not occur on 
circular orbits and a more detailed calculation is required. 



Gas-drag induced capture and coagulation of 
planetesimals: Two unbound planetesimals may dis- 
sipate some of their kinetic energy during an encounter 
(due to gas-drag), and may become bound to form a 
transitional binary. Such binaries could then continue to 
inspiral and finally merge due to gas drag (as discussed 
in the previous section). This capture-coalescence pro- 
cess, which could play a role in the build up and coag- 
ulation of planetesimals will be discussed in detail else- 
where (Murray-Clay and Perets, in preparation; see also 
the settling regime discussed by Ormel & Klahr, 2010). 

Planetesimal erosion: As shown above, in some 
regimes the WISH radius is smaller than the big plan- 
etesimal size. This would suggest that under these condi- 
tions wind-shearing may blow away loose parts from the 
surface of single planetesimals, if they are weakly bound 
to the planetesimal (e.g. pieces from aggregates held to- 
gether only by gravity could be blown from the surface of 
the main component of the planetesimals). Indeed, such 
erosion of planet esimals was experime ntally observed for 
dust aggregates (IPa raskov et al.|[2006[ ). 

Post impact evolution of planetesimals: Follow- 
ing the collisions of two planetesimals some of their ma- 
terial may be ejected from the surface. A large frac- 
tion of this material is still gravitationally bound to the 
system. However, small size particles gravitationally 
bound to larger planetesimals could be blown away by 
the wind, if they are ejected beyond the WISH radius. 
The post-impact evolution of these particles could there- 
fore be qualitatively different than the corresponding 
non-gaseous collisional evolution, prohibiting the small- 
est impact debris particles from ever accreting to the 
main bodies of the planetesimals. Nevertheless, wind 
may also induce re- accretion of ejecta material in some 
cases (Teiser & Wurm 2009), depending on the ejecta 
trajectory. In addition, the short merger timescales we 
find for binary planetesimals, suggest that any (WISH- 
stable) bound debris around the main collision remnant 
would inspiral and accrete to the main body. 

Planetesimal encounters and collisions in a 
gaseous environment: The collision rates between 
planetesimals vary for different regimes of velocity and 
encounter distances scales, where one of the most im- 
portant scales of the problem is the Hill radius (e.g., 
lGoldreichet"aL|[200l . The WISH radius provides an 
additional important parameter for planetesimal encoun- 
ters, which has to be taken into account in order to de- 
termine the outcome of plan etesimal encounters in gas. 
Recently (and independently) lOrmel fc Klahrl ()201ClD dis- 
cussed planetesimal encounters in a gaseous environment 
for some specific encounter regimes, and provided de- 
tailed calculations for these regimes. Our study suggests 
an additional and complementary understanding of these 
issues. 

5. SUMMARY 

In this study we explored the differential gas-drag 
acceleration between different size planetesimals in a 
gaseous environment. We defined the wind-shearing ra- 
dius as the distance at which the differential accelera- 
tion between two close-by planetesimals is comparable 
to their mutual gravitational pull. Planetesimal inter- 
actions close to or beyond this limit would be strongly 
affected by wind shearing. The wind-shearing radius has 
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important implications for the existence and survival of 
binaries. We find that binary planetesimals cannot form 
or survive with separations beyond this scale, even if this 
separation is smaller than the Hill radius, as they would 
be destabilized and sheared apart by the head wind. 
WISH-stable binary planetesimal are also affected by gas 



drag, and can inspiral and coalesce in times shorter than 
the lifetime of the gaseous disk. The wind-shearing ra- 
dius may have important implications for planetesimal 
evolution, in particular planetesimal erosion, post impact 
evolution of planetesimals, and planetesimal encounters 
and coagulation. These effects merit further investiga- 
tion. 



APPENDIX 

WIND-SHEARING DISRUPTION OF BINARY PLANETESIMALS : PARAMETER DEPENDENCE 

In the following we provide a more detailed discussion on the dpendence of the WISH stability criterion for binary 
planetesimals on their environment and properties. 

Gas density and relative velocity 

The gaseous environment and the velocities of planetesimals with respect to the gas likely change with time in a 
given disk, and vary across planetary systems. Therefore, the WISH radius and the stability of binary planetesimals 
and satellites is time and system dependent. These issues are illustrated in Figure [5l which shows how the WISH 
radius varies as a function of the relative velocity and gas density. The figure shows the WISH stability radius of an 
Tb = 10 km planetesimal at 1 AU, for various choices of the relative planetesimal-gas velocity (upper panel) and the gas 
density (lower panel) . Note that for small planetesimals in the Stokes regime the WISH radius becomes independent of 
the gas-density, as can be seen in Eq. ([T^ . and the various lines in Figure [S] converge. We do not show the dependence 

of Rws on temperature since it is weak (at most Rws oc Cs oc T^i/* ; gee Eq. [T^ . 

As might be expected, the general trend of the WISH radius is to be smaller for higher gas densities and/or higher 
planetesimal velocities relative to the gas, i.e. WISH becomes more pronounced with stronger gas drag. Also, as men- 
tioned before, very large comparable size planetesimals are hardly affected by gas-drag and the WISH radius becomes 
larger than their Hill radius. Taken together, the WISH will be more important for smaller binary planetesimals, 
during earlier stages of their growth/evolution, in a more gas-rich environment; WISH gradually becomes negligible 
at later stages. 
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Figure 5. Dependence of the the binary WISH stability radius on the gas-density (top) and the relative velocity between the planetesimals 
and the gas (bottom). These figures show the dependence for a specific choice of the big planetesimal, rj, = 10 km at 1 AU in our fiducial 
disk. The density and velocity for each line in the upper and lower panels, respectively, are S = /s X 2 X 10^ g cm"'^ and v^ei = fvX CSc^/i^x. 
The dashed lines in both panels correspond to the same parameters used for the respective lines for the 10 km size big planetesimals shown 
in Figs. [T]and[2] Shaded regions correspond to separations beyond the Hill radius (upper region) or below the physical radius of the big 
planetesimals (lower region). 



